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Abstract 

We show how to build up a systematic bosonic field-theory for a general reaction-diffusion pro- 
cess involving hard-core particles in arbitrary dimension. We criticize a recent approach pro- 
posed by Park, Kim and Park (Phys. Rev. E 62 (2000)). As a testbench for our method, we show 
how to recover the equivalence between asymmetric diffusion of excluding particles and the 
noisy Burgers equation. 
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1 Introduction 

A general formalism to describe nonequilibrium statistical physics is still lack- 
ing, but recent progress has been achieved on the important issue of the description 
of collective phenomena in such systems (phase transitions in nonequilibrium steady- 
states [p], emergence of long-range correlations in driven systems [ffl], cellular au- 
tomata [3], self-organized criticality [||]). Indeed steps forward have been made in 
the last two decades in that direction by exploiting the formal analogy of their field- 
theoretic formulation with static critical phenomena. Of course, numerical approaches 
and exact solutions have played an equally important role, but this communication will 
deliberately ignore those aspects. The beauty of the field-theoretic formulation -when 
available- is that, combined with a renormalization group analysis, it provides the the- 
orist with a systematic analytic tool for the calculation of physical observables in the 
scaling regimes of interest. However, a vast category of reaction-diffusion processes, 
cellular automata, driven lattice gases and other related stochastic models fail to be ex- 
actly mapped onto continuous field theories. This difficulty is due to the fact that the 
integer degrees of freedom (often particle numbers) usually have an integer upper limit 
which prevents from exploiting the familiar mapping first introduced by Doi [Q] and re- 
cently revived by Cardy [||| . For instance, in a reaction-diffusion process with mutually 
excluding particles, local particle numbers are or 1, but one might imagine other pro- 
cesses with other similar constraints (e.g. the contact threshold transfer process defined 
by Rossi et al. ^ in which a given site is occupied by no more than two particles). In 
some instances, the exclusion constraint can be phenomenologically accounted for, as 
done e.g. by Zia and Schmittmann [^] when they build up effective Langevin equation 
to describe the dynamics of driven-diffusive systems. Very often the evolution operator 
can be exactly mapped on a quantum spin-chain. But it is only in one space dimension 
that (in favorable cases) one can exploit the toolbox of integrable systems [n0| (local 



particle numbers usually restricted to and 1, little being known [11 1 on spin-1 (and 
higher) chains for stochastic problems). 



A coveted goal is therefore to be able to build up a systematic and exact field- 
theoretic path-integral formulation that can, at least formally, account for limitations in 
local particle numbers. 

Several attempts have been made to incorporate the exclusion constraint into a 
field-theory of standard type. As already mentioned, one may achieve this goal by phe- 
nomenologically exploiting the physical knowledge of the system (Zia and Schmittmann 
or Cardy [p"2[|), but this is by no means exact, and can in general only be implemented 
safely in high space-dimension. Directly in dimension 1, other approaches exist, such 
as that proposed by Cardy Jl3| ] and Brunei al. [ |l4] ] or by Mobilia and Bares [Q. In the 
former the authors succeeded in constructing a field-theory in a systematic way, but of 
fermionic type, which unfortunately proved, from a technical view, rather difficult to 
analyze, while the latter, though efficient in the pair annihilation reaction, seems diffi- 
cult to extend to other processes. 

In what follows we shall present the derivation of the path-integral formulation for 
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systems of mutually excluding particles of a single species. Where appropriate we will 
also indicate how to extend the theory to cope with several species or other constraints 
on local particle numbers. We shall illustrate how the mapping works on the case of 
asymmetric diffusion, and how one can recover the equivalence with the noisy Burgers 
equation. 

2 Hard-core particles using a bosonic formulation 

2.1 Master equation and bosonic formalism 

The evolution of a configuration n = {rij} of local particle numbers is encoded in 
a master equation for the probability P(n, t) to observe configuration n at time t. The 
master equation for P(n, t) is equivalent to an evolution equation for the state vector 
\fy{t)) = J2 n P( n i t)\ n )i which we write in the form 

d t \*) = -£|*> (i) 

The operator H, which acts on the space spanned by the configuration vectors |n) is 
usually easily expressed in terms of bosonic creation and annihilation operators a\ , a% 
= Sij, [ai,a,j] = 0). This is true for reaction-diffusion processes involving 
bosonic particles, that is, without exclusion, for which bosonic operators are particu- 
larly well-suited, but this is also true when particles exclude each other. Nevertheless 
in the latter situation one should not expect that H will be a polynomial in terms of the 
a], while it indeed is in the former. 

We confine the subsequent analysis to the dynamics of hard-core (i.e. mutually 
excluding particles) undergoing diffusion and reaction processes. 

2.2 Evolution operators for some elementary processes 

In order to deal with the exclusion constraint we introduce the operator Sa m defined 

by 

$n, m \m') = 5 m > im \m') (2) 

where \m) denotes a single site state. We now give two examples. The evolution opera- 
tor for one-dimensional asymmetric diffusion in the presence of hard-core interactions 
reads 

H m = ^2 [( D + l^ 1 _ a »°]+i) < W < Wi,o + ( D ~ |)(! - a t Wi)^,o<5ft 1+1 ,i 

i 

(3) 

where D + | (resp. D — |) is the hopping rate to the right (resp. to the left). 
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For the simple A + A — > annihilation reaction of nearest neighbor particles, one 
finds 



H k = [(! - 4 a i+i) < W<Wi, 



(4) 



Extension to two-species annihilation is straighforward since bosonic operators per- 
taining to distinct species commute. 



2.3 Passing to a coherent state representation 

In order to pass to a path-integral formulation, it is sufficient to follow the steps de- 
scribed in [ |l6| l (for a thorough and pedagogical introduction we refer the reader to the 
review by Mattis and Glasser [|r7[|). The result of those steps can be summarized as fol- 
lows. There exists an action S[cj), <fi] such that physical observables can be expressed as 
path-integrals over the complex fields (j>i(t), 4>i(t) of functions of those fields, weighted 
by exp(— S). We denote by [0, tf] the time interval over which the process is studied. 
The action S has the form 

S=-J2<f>i(tf)+ r dt (Y.MWth+H&t]) (5) 

where 

TJU M WjM (to 

in which the notation \<j>) — ®i\fa{t)) denotes the tensor product of the coherent states 
associated to each creator and annihilator a\ , di with eigenvalue fa (t) (and fa (t) de- 
notes the complex conjugate of fa(t)). In order to evaluate the quantity ^j^j,^ , one 

normal orders H and then simply replaces the ct^'s by fa(t) and the aj's by fa(t). 
Hence the only difficulty is to be able to normal order such an operator as 5a . m (pos- 
sibly multiplied by a's or a''s). This can be done rather easily in several ways. One of 
them is simply to normal order e lun then make use of the integral representation 

x du 



X. — I ^ U Jujh-m) 
fJ n,m / c 

J-7T 27T 



Normal ordering e zun is done by expanding the exponential and looking at each term 
in the series. It is amusing to note that 



E^E 

1=0 ' j=l 



Sj^a^a-' (7) 



where the coefficients Sj.e are the Stirling numbers of the second kind (sjj = 4^- (e x - 
iy is the number of ways a set with I elements can be partitioned into j disjoint, 

z—0 
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non-empty subsets). Similarly, one finds that 



71 ! 
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Once in their normal-ordered form, the operators a and in Eqs .([7]) and (Q) can be 
replaced by their coherent state eigenvalues. For instance, 

+ 00 I ■ \t & 



so that 

(rh\X* n \rh\ = /" 

„ 2tt 

Similar manipulations lead to the dictionary 

(4>\a)5f r 



(0nM= r ^e^- 1 )=e-^ (10) 
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to 



to! 

MaS^M = 4> , 1 1M Wr-V^ (=0ifm = 0) (11) 
(to — lj! 

In a reaction-diffusion process with hard-core particles, only the formulas for to = 
or 1 will be needed. Formulas with higher powers of the a's or a^'s are derived in a 
similar way. 

For example, the action for the annihilation process described by the evolution 
operator of Eq. (Eh reads 



S = k / dtE ~ l)&0t+i e 



(12) 



Extensions to several species poses no new problem. A final remark : when studying a 
process that conserves the parity of the instantaneous total number of particles N(t) in 
the system, the following quantity is conserved 

((-1)^*)) = ( e i7rAr «) = (p|e" rS * fl *|1<) (13) 

where (p| = (0|e^-- ^ ai denotes the projection state [|6|, which is also a coherent state 
with eigenvalue 1, so that, using that (l|e I7mi \<f>i) — e*^ 6 " -1 ), we find 

((-l)^) = ( e - 2 S«*i(*)) =cs t (14) 

This is a way of characterizing parity by means of a well-defined observable. We refer 
the reader to Deloubriere and Hilhorst [ |l8| ] for further comments in the context of the 
pair annihilation reaction. 
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2.4 The Park, Kim and Park approach 

It is well-known [ |l9| ] that for reaction-diffusion processes involving bosonic parti- 
cles (which is not the case in it is possible to write a partial differential equation 
for some continuous random variable pi (the mapping fails for hard-core particles). 
The first moment of pi equals the average local particle number (rii) (which makes 
it tempting to identify pi with a fluctuating density). However higher moments of pi 
do not coincide with those of m (though they can be related). This partial differential 
equation takes a Fokker-Planck form (i.e. is of order 2) only when the microscopic re- 
actions involve at most two particles. 

We now refer to the article [[[J in which the authors have presented an alternative 
route to derive a path-integral formulation for the dynamics of hard-core particle sys- 
tems. Particle numbers are discrete variables, so that there is no Fokker-Planck equation 
for them. In their Eq. (7) they write a Fokker-Planck equation for some continuous ran- 
dom variable pi, which, as we have said, is not correct without further approximation. 
This error is not related to the necessity of implementing the hard-core constraint or 
not. Such a Fokker-Planck equation simply does not exist. 

References and further comments can be found in the book by Gardiner [|l9|] or in 
Deloubriere and Hilhorst @. 

3 An example : Asymmetric diffusion of hard-core par- 
ticles and the noisy Burgers equation 

As an example we recover the noisy Burgers equation by going to the continuous 
limit in the asymmetric diffusion of a system of hard-core particles. While this equiv- 
alence is certainly not new f?~2"| l, we use it as a testbench for the method. It 
should be mentioned that this derivation of the noisy Burgers equation is the first one 
that starts from an exact mapping. 

3.1 Action 

Again, for notational simplicity, we restrict the analysis to one space dimension. 
We consider diffusion with a hopping rate to the right D + % and a hopping rate to the 
left D — |. The evolution operator for such asymmetric diffusion in the presence of 
hard-core interactions is that of Eq. ([}]). Using the dictionary Eq. ( |Tl| ) we find that the 
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corresponding action reads 



^ +1 )^e-^^-^+ 1 ^+ 1 (15) 
^xj&e-*'*-*'- 1 *- 1 

No approximation was made and the action Eq. ( [l5] ) is fully exact. Without the ex- 
ponential factors Eq. ( |l5| ) would yield the usual action of asymmetric diffusion for 
bosonic particles. Here, owing to the presence of the nonlinear interaction terms (the 
exponentials), there is no Galilean transformation that eliminates the drift-dependent 
terms. 

3.2 Recovering the noisy Burgers equation in the continuous limit 

In this paragraph, we shall show that expanding naively the action Eq. ( |l5| ) leads 
to the noisy Burgers equation for the density fluctuations, as it should. We perform the 
change of fields [ p3| fa = (p + fa)e~^\ fa = e^\ Now we expand the action in 
powers of the new fields fa and fa (the latter represents a fluctuation of the density 
with respect to its average value p). We also assume that the fields have slow space 
variations, and take the limit of a continuous space. The resulting action reads 

S = J dtdx [{b{d t + vd x - Ddl)ip - gi{d x ^f - g 2 ^d x ^ + ...] (16) 

where the constants gi , g2 are positive functions of the microscopic details of the model 
(diffusion constant, drift velocity, lattice spacing, average density). The dots stand for 
higher polynomial or higher derivative terms. The action Eq. (|lq ) for the fields ?/;, tp 
is equivalent to ip(x, t) satisfying the noisy Burgers equation []22[]. Hence, up to terms 
which are irrelevant in the scaling limit (see Janssen and Schmittmann p2| ] for a renor- 
malization group analysis of the action Eq. (p"oj)), we have recovered the equivalence 
between asymmetric diffusion of hard-core particles and the noisy Burgers equation. 
As a final remark, we would like to emphasize that we have not resorted to any phe- 
nomenological arguments, and that our derivation is completely systematic - the first 
one of this sort. This feature is particularly encouraging since we have in mind the 
application of the formalism to other less studied processes. 

4 Conclusions 

We have shown how to build up a field-theoretic formalism that takes into account 
in a systematic fashion the effect of exclusion in d-dimensional reaction-diffusion pro- 
cesses. We have exemplified the formalism on the case of asymmetric diffusion, thus 
recovering the noisy Burgers equation. We now have a tool to take up any reaction- 
diffusion process in which one, or all species, diffuse with a drift, such as the A + B — > 



sr. 



<P + -)(4* 



(D--)(fa 
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reaction, for which it is conjectured that exclusion changes the universality class of 
the scaling behavior. The A + B — > reaction-diffusion process with drift is certainly 
the first system that should be looked at using the present approach. However, great 
care must be paid to naive expansions of exponential interaction terms, and the feasi- 
bility of such a procedure must be investigated in each particular case. 

The method presented here opens the door to the study of reaction-diffusion pro- 
cesses in which exclusion is conjectured to play a crucial role, such as in the N species 
branching annihilating random walks recently described by Kwon et al. [[24]]. 

Acknowledgments : The author would like to thank Henk Hilhorst, Uwe Tauber and 
Martin Howard without whom this approach would not have been elaborated, and 
would like to thank Gunter Schiitz for communicating [ |To| ] prior to publication. 
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